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Abstract 

A study at both the molecular and extended solid level in the framework DFT is carried out for KM[Cr(CN)e] ( M = V, Ni). 
From molecular calculations, the exchange parameters J are obtained, pointing to the expected magnetic ground states, i.e., 
antiferromagnetic for M = V with J = -296.5 cm -1 and ferromagnetic for M = Ni with J = +40.5 cm - . From solid state 
computations the same ground states and J magnitudes are confirmed from energy differences. Furthermore an analysis of the site 
projected density of states and of the chemical bonding is developed in which the cyanide ion linkage is analyzed addressing some 
isomerism aspects. 
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1. Introduction 

Prussian blue KFc 111 [Fc 11 (CN)g] , is a historically known 
chemical compound pQ. The presence of two types of 
cationic transition metals (TM) as shown in fig. [Hallows to 
extend the family to several systems by substitutions not 
only of TM but of K as well by other monovalent alcalinc 
cations such as Cs [2] ■ Although the crystal structure is face 
centered cubic with high symmetry (space group F-43m, 
N°216), the elucidation of the Prussian blue structure was 
the object of intense research work due to the crystallo- 
graphic problems [3]. Nowadays, the study of magnetism in 
such systems is a rich field of investigation and many mag- 
netic materials have been synthesized and characterized. 
Furthermore, due to the high symmetry of the structure, 
molecular field theory can be easily applied to these types 
of systems [J . Therefore, Prussian blue analogues are the 
obvious choice for the design of molecular magnets, and 
many fascinating and exciting phenomena have been re- 
ported in the literature. Examples of such studies include 
high ordering-temperature magnets |5|6|7j . photoinduced 
magnetism [5] , electrochemically tunable magnets [Hj , and 
magnets that experience two compensation temperatures 
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[ID] . Most of these examples have been designed taking 
advantage of the properties of this family of compounds. 
The solids of Prussian blue analogues are readily prepared 
from cyanometallate blocks and transition metal cations. 
The first quantum theoretical studies were done on the 
molecular systems using the extended-Hiickel method for 
dinuclcar models |5|llj . They applied the Kahn-Briat 
model [12] and attributed the exchange interaction to the 
overlap between magnetic orbitals through the tt system. 
Many other studies on the exchange interactions were 
carried out to evaluate the magnitude of J coupling pa- 
rameter allowing an assignment of the type of magnetic 
order from molecular calculations [13|14j . For the solid 
state periodic structures, the work of Eyert et al. [15] on 
two well-characterized Prussian blue analogues of formula 
CsM[Cr(CN) 6 ] (M = Mn, Ni) can be noted within density 
functional theory (DFT) [l"6] . 

In this work, we carry out a study at both the molecular 
and extended solid levels within DFT for KM[Cr(CN) 6 ] (M 
= V, Ni). Exchange parameters J, are obtained from molec- 
ular calculations and confronted with J values from en- 
ergy differences within solid state computations assuming 
different magnetic configurations, namely non-magnetic 
(NM), ferromagnetic (FM) and antiferromagnetic (AFM) 
configurations. Special attention will be devoted to an 
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analysis of the site projected density of states and a study 
to the chemical bonding is developed in which the cyanide 
ion linkage is analysed addressing isomerism aspects. 

This double approach at the molecular and solid state 
levels has already allowed us to address the low spin 
and high spin properties of the spin cross-over complex 
[Fe(btz) 2 (NCS) 2 ] PI]. 

2. Magnetic orders 

Depending upon the response of a material to an applied 
field, we can describe it as either a diamagnetic or para- 
magnetic sample. Paramagnetic materials may be further 
classificld as fcrromagnets, fcrrimagncts and antifcrromag- 
nets. Table I sketches out the ferromagnetic superexchange 
interactions through the cyanide bridge between d orbitals 
of Cr and M with environment. This diagram involves 
a difference of electron occupation of e g which can lead to 
ferromagnetic or antiferromagnetic interactions. When the 
interaction is between t 2ff - e s which is the case of Cr-Ni, 
this leads to a ferromagnetic interaction due to the elec- 
tronic structure of the cyanide bridge and special role of 
nitrogen, while t2 S Cr - t2 9 V leads to antiferromagnetic in- 
teractions [15] , 

Fcrromagnets, behave as paramagnets at high tempera- 
tures, but below a critical temperature Tc (Curie temper- 
ature), they show magnetic ordering. In a mean field ap- 
proximation, it can be demonstrated that: 

\/ ZmZk[>\J\^JxSm{Sm + 1)Sm'(Sm' + 1) HS 
TC = 3k~ B (1) 

where Sm and S M are the local spins on centers M and 
M , Zm and Z M the number of nearest neighbors of each 
type of metal atom and x corresponds to the stoichiomctry 
of the M'[M(CN) 6 ]:r compound. J is the exchange coupling 
between the two TM. It indicates the type of communi- 
cation between the two spin centers. A positive value of J 
represents a ferromagnetic coupling, while negative J de- 
scribes an antiferromagnetic interaction. The J values are 
obtained by using equation: 

E HS -E LS = -(2S 1 S 2 +S 2 )J (2) 

where Si and S2 are the total spins of the paramagnetic 
centers and Si > S2 has been assumed for heterodinuclear 
complexes. E#s is the energy of highest spin state and E^g 
that of the sate having the two metals with opposite spins. 
Over the last decades, the Prussian blue analogues 
A/M m [M(CN)g] (H2O), have received attention as molecule- 
based magnetic materials. Effective exchange interactions 
between M and M via the CN~ anion (note that the neg- 
ative charge is carried by carbon strong leading to high 
Tp values. Verdaguer and coworkers have synthesized 
compounds with Tc value of 90, 240 and 315 K, with M 
= Ni(II) and M = Cr(III) [51X5] . M' = Cr(II) and M = 
Cr(III) [20], and M' = V(II) and V(III), M = Cr(III) [Jj 
respectively. 



3. Computational Details 

3.1. Calculations at the molecular level 

Although the use of Hartree-Fock (HF) approach has 
been shown to provide a good description of the molecular 
orbital and chemical bonding properties mainly in organic 
chemistry, it becomes well established that calling for the 
density functional theory (DFT) framework |16| brings far 
more accurate results regarding the energetics and related 
properties. This is because the compulsory exchange and 
correlation (XC) effects are equally treated, albeit at a local 
level, within DFT while only exchange is well accounted for 
in HF although in abetter way (exact exchange [21]) than in 
DFT. Taking the best out of each one of the two approaches 
led to improvements in ab initio molecular calculations with 
the so called " hybrid functionals" . They consist of mixing 
exact HF exchange, ex. following Bccke [21] and DFT based 
correlation, ex. following Lee, Yang and Parr, i.e., the so- 
called LYP correlation [22 j . with proportions that help to 
reproduce molecular properties of several systems. The cal- 
culations have been performed using GaussianOS package 
[23] . Owing to the complexity of the electronic structure of 
the studied compounds, we have employed the JAGUAR 
code [24| to generate an appropriate initial wavefunction set 
since it provides a better control of the local spin and multi- 
plicities of the atoms. To maintain the electronic structure 
during the SCF procedure the quadratic convergence op- 
tion must be used. In all calculations, the hybrid functional 
B3LYP and the LanL2DZ basis sets were used. The latter 
includes double-^ with the Los Alamos effective core poten- 
tial for Cr, V, Ni and the Dunning-Huzinaga all-electron 
double-^ basis set with polarization functions for the H, 
C, and N atoms (£ is the exponent in the Gaussian type 
orbitals GTO) [25|26j . The following bond lengths used 
for the calculations, were for the KV[Cr(CN)g] complex: 
dc r -c = 2.01 A, dv_7v = 2.01A and d C - N = 1-24 A. For 
the second complex the distances were : 2.03, 2.04 and 1.24 
for dcr-cs d-v-N and dc-jv respectively. 

3.2. Solid state computations 

In as far as no structural determinations were available 
for both systems, trends in cell volumes and ground state 
crystal structure were desirable in the first place. The equi- 
librium structures were obtained from a full geometry re- 
laxation from CsNi[Cr(CN)g] structural setup [15] and us- 
ing a pseudo-potential approach within the VASP package 
[27] . Ultra-soft Vanderbilt pseudo-potentials (US-PP) built 
within the local density approximation (LDA) scheme [28] 
were used. The Brillouin-zone integrals were approximated 
using a special fc-point sampling. These values were then 
used in the all electron augmented spherical wave method 
(ASW) [25]. Besides its use of DFT, this method is based 
on the atomic sphere approximation (ASA) , a special form 
of muffin-tin approximation which consists in dividing the 



cell volume into atomic spheres whose total volume is equal 
to the cell volume. The calculation being carried out in 
the atomic spheres space, empty spheres (pseudo-atoms) 
need to be introduced -without symmetry breaking- in open 
(low compactness) structures such as that of the systems 
studied here. Further, empty spheres allow for the iono- 
covalent characteristics of the system to be accounted for 
by receiving charges from actual atomic species. The use of 
this method in molecular systems was formerly validated 
by Eyert et al. [15] . The Brillouin-zone sampling was done 
using an increased number of k points ranging from 126 to 
840 points within the irreductible wedge. In this way we 
were able to ensure convergence of the results with respect 
to the fineness of the k-space grid. 

In the context of this study, ASW-LDA was used in order 
to obtain a description of the electronic structure with the 
partial, site projected, density of states DOS. Further the 
chemical interaction is assessed by using the ECOV (cova- 
lence bond energy) criterion, which allows to get two-body 
chemical bond characteristics [30]. Negative, positive and 
zero ECOV magnitudes point to bonding, anti-bonding and 
non-bonding interactions respectively. In order to enable 
for a better understanding of J, the electronic properties as 
well as the DOS and the ECOV, we start out by discussing 
the results for both compounds. The strategy of the ASW 
calculations was to do three types of computations for the 
two complexes : Firstly a non-magnetic NM calculation 
(degenerate spins) , then detecting an instability of the NM 
configuration magnetic calculation, a Ferromagnetic calcu- 
lation is carried out. Finally, the system is enforced to be 
anti-ferromagnetic through supercell calculations in order 
to establish the ground state expected for KV[Cr(CN)g]. 

4. Results and discussion 

4.1. Calculation of J 

The three transition metals, chromium, vanadium and 
nickel, ionize as Cr /7/ (d 3 ), V H (d 3 ) and Ni /7 (d 8 ) within 
the complexes KM[Cr(CN)g] under consideration. These 
formal degrees of ionization are fed into JAGUAR input in 
order to calculate the charge densities of the different mag- 
netic alignment configurations, i.e., FM and AFM. Then 
with unrestricted HF (accounting for spins) calculations 
the energy differences (cqn. [2]) allow obtaining J. For M = 
V system, a very strong antiferromagnetic coupling with a 
large J = -296 cm -1 value was found. This corresponds to 
the scheme of V /7 -t2 9 < — > Cr //7 -t2 9 antiparallel coupling 
with Oh environment, in both e g manifolds are empty (see 
lower part of Table 1). On the contrary a ferromagnetic 
coupling with J = + 40.5 enr 1 value is found to agree with 
the Cr //7 -Ni // parallel coupling of spins: tf> g < — ► ^2g e g- F° r 
these two complexes, Ruiz et al [13] have obtained the val- 
ues of J = -241 cm -1 and +22.4 cm -1 . From this outcome 
our values follow the same trends but the differences can be 
due to the variation of the distances and the basis set they 



use, i.e., TZV ( triple £). Tq calculated via eqn. [TJ gives a 
value of 3171 K (for 315 K) for the KV[Cr(CN) 6 ], and 274 
K for KNi[Cr(CN) 6 ] (for 90 K) [13]. This overestimation of 
the Tc values is a well-known problem. The electron con- 
figuration of M and M' affects the nature and magnitude of 
the coupling in a predictable way, t^-t^ pairs giving strong 
antiferromagnetic interaction (case of Cr-V) while t 2 ™ -t\ 
pairs yielding ferromagnetic coupling (case of Cr-Ni) (See 
Tabic 1). This result is in agreement with those obtained 
by Ruiz [13] who have proposed these explanations based 
on the calculations with several transition metals. Table 2 
regroups the atomic spin densities on Cr, V and Ni calcu- 
lated with the hybrid functional B3LYP and the LanL2DZ 
basis sets. The density of state on Cr varies slightly from 
ferromagnetic to the antiferromagnetic state. 
In addition to the calculation of J, the infra-red (IR) and 
Raman theoretical frequency spectra were computed by us- 
ing also the hybrid functional B3LYP and the LanL2DZ 
basis sets cf. fig. [2] The most important peaks are those 
at v =189.017 enr 1 assigned to the V-N stretching, Cr-C 
elongation appears at 255.447 enr 1 , the high intensity at 
1835.03 cm" 1 is that of elongation of the C = N between 
the two metals. The latter is also found in Raman spec- 
trum. The elongation of the other C=N appear at 2051.27 
cm" 1 as an intense peak. 

The frequency of the electromagnetic wave which induces 
the vibration of elongation is given by the relation 

where k is the constant of bonding strength (considered 
here as a spring) , proportional to the binding energy and /i 
the reduced mass of the two atoms connected by this bond. 
Feeding the v values in eqn.[3]for Cr-C and V-N provides a 
ratio of force constants k, as k( y r ~jv ) = 1-^2 which means 
that the strength of Cr-C bond is 1.62 larger than that of V- 
N bond. This will be confirmed in the study of the chemical 
bond in the solid state. 



4.2. ASW calculations 

KV[Cr(CN) 6 ] and KNi[Cr(CN) 6 ] both crystallize in a 
face-centered cubic lattice with space group F43m [6]. In 
this space group K, M and Cr are found at fixed Wyck- 
off sites, i.e. respectively at (4c/4d), (4a) and (4b) while 
the (24f) position at x,0,0 allows for variable x value. The 
lattice relaxation with VASP US-PP calculations keep the 
fee symmetry and give the following values for (24f) posi- 
tions: X p /Wl = 0.309/ 0.306 and = 0.192/0.188. A 
good agreement with literature [5] is found for the C and 
N x positions for the M = Ni, Cs-based system. The re- 
laxed lattice constants are &KV[Cr(CN) 6 \ = 10.52 A and 
a -KNi[Cr(CN) 6 ] = 10.64 A which come in close magnitudes 
to other similar systems |15] . 
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4.2.1. Nonmagnetic calculations for KV[Cr(CN)§] 

At self-consistent convergence of the energy and of 
charges, the charge transfer is found from TM towards C 
and N as well as to the ES with a larger magnitude. In this 
hypothetical magnetic configuration the system is found to 
behave as a metal due to the large density of states at the 
Fermi level. This is exhibited in the site projected partial 
densities of state (PDOS) in fig. [4ji illustrating the results 
for KV[Cr(CN)g] member. All energies are relative to the 
Fermi level . In the closely octahedral (O^) crystal field 
of CN ligand surrounding the TM ions, V and Cr, the d 
orbitals split into two main manifolds which are labeled as 
t2 g (three-fold) and e g (two-fold). In the PDOS plots three 
dominant peaks are identified for Cr at -4 eV for the states 
which mix with CN ligands, at Ej? for non-bonding t2 S 
orbitals and at 3.5 eV, in the conduction band, for empty 
e g . Vanadium states are characterized similarly with DOS 
peaks at -6 eV, at Ep for t2 g and at 2.5 eV for empty e g . 
Carbon and nitrogen DOS in the lower energy part of the 
valence band, i.e. in the [-6,-2 eV] energy range are very 
similar pointing to the covalcnt C...N bonding, in fact the 
cyanide is built through strong interactions (one a and two 
it) between carbon and nitrogen orbitals. 

The very high density of states at the Fermi level char- 
acterizing V and Cr is a sign of instability for the system 
in such a NM configuration. This is analyzed through the 
Stoner mean field theory of band ferromagnetism |31j . If 
the energy of the system is counted from a NM configura- 
tion, it can be expressed as follows: 



In this equation, m is the magnetic moment, I is the Stoner 
exchange-correlation integral which is an atomic quantity 
that can be derived from spin polarized calculations, tief 
represents the PDOS value of the element under considera- 
tion at the Fermi level in the NM state. Between the square 
brackets in eqn.[4j the product I.yief can be considered as 
criterion of instability of the system: if it is >1, then the 
energy counted from NM configuration is lowered because 
1-I.iiEF < 0, hence the system is stabilized by magnetic ex- 
change. The product I.n^ in cqn. 3]is then considered as 
a criterion for the stability of the spin system. For V, I.tief 
= 1.831 and in the case of Cr it is equal to 1.900. From this, 
the system is unstable in a nonmagnetic configuration. 

We now turn to addressing the chemical bonding through 
the ECOV criterion introduced above. In fig.[5jD, along the 
y-axis, unitless positive, negative and zero ECOV magni- 
tudes point to bonding, anti-bonding and non-bonding in- 
teractions respectively. The bonding between the TM ions 
and C/N is mainly ensured by the itinerant-like electrons 
in the lower part of the valence band, i.e., in the energy 
range [-4,-2 eV] where one observes that the Cr-C bonding 
is twice larger than the V-N one. This confirms the finding 
of the ratio k( V;J^y ) from the vibrational spectra analysis 
of the isolated molecules providing the same trend without 
having the same value. At the Fermi level there are mainly 



t2 9 -like orbitals which show antibonding character; but this 
is not the major issue of the bonding within the system 
as these orbitals are mainly engaged in the magnetic in- 
stability as discussed above. In the conduction band above 
F>f the antibionding counterparts of the Cr-C and V-N are 
found in the energy range [2, 6 eV]. 

Furthermore it can be suggested from the analyses above 
from the antibonding interaction between V and C as well 
Cr and N that V must be bonded to N and not C, con- 
trary to Cr which must be bonded to C. In order to pro- 
vide rationalize this observation, an additional calculation 
was carried out with the hypothesis of Cr-N and V-C di- 
rect bonds. The resulting total energy is found to desta- 
bilize the system by 1.32 eV. Thus the system is stable in 
the experimentally proposed configuration. This result is 
further confirmed using pseudo-potential calculations with 
the VASP code, where the destabilization amounts to 1.28 
eV. Such energy criteria are features that can be used as a 
predictive tool for unknown systems. 

4.2.2. Magnetic calculations for KV[Cr(CN) 6 J 

From spin occupation for each atom species at self con- 
vergence, the total magnetization is close to 6 for the unit 
cell. This agrees with Hunds rules' spin-only magnetic mo- 
ments of 3 for both V and Cr. From the partial DOS cf. fig. 
2k, the spin polarization clearly affects the two TM Cr and 
V. Also, we have a very low DOS at the Ey pointing to the 
top of valence band. Total energy favors the ferromagnetic 
spin state state by 1.15 eV versus the non magnetic config- 
uration. From Table 3, all other atomic species, including 
the empty spheres exhibit negligible moments. 

4.2.3. Antiferromagnetic calculations 

Using a double cell with anti-parallel directions, the to- 
tal magnetization is equal to zero and the system becomes 
an insulator with a gap opening of ~ leV cf. fig. [4)d. The 
spin magnetic moments for the ferro- FM and antiferro- 
magnetic AFM configurations are given in table 3. Their 
magnitudes show a small lowering with respect to the FM 
case; this is probably due to the symmetry breaking. To- 
tal energy favors the AFM state by 0.154 eV versus FM. 
This points to the antiferromagnetic ground state for the 
KV[Cr(CN) 6 ] complex. 

The calculated J value in the solid state from the same 
equation is equal to -418 cm -1 which is found 41% larger 
than the molecular calculations (-296.5 cm -1 ). Such differ- 
ences could be expected due to the different account of the 
chemical system as an isolated molecule versus an extended 
solid; nevertheless the trends are the same. 

4.2.4. Calculations for KM[Cr(CN) 6 ] 

Following the same calculation protocol, the nonmag- 
netic configuration is found unstable to the ferromagnetic 
one by 1.6 eV. From the magnetic calculations, the total 
magnetic moment of the cell is equal to 5 in agreement 
with Hund rule Ni 2+ and Cr 3+ . The atomic magnetic mo- 
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merits are given in the table 5. Contrary to the case of 
KV[Cr(CN)@] a hypothetic AFM configuration is found less 
stable than FM thus confirming the magnetic ground state 
for the system (see table 4) . 

The energies of the two complexes in the three configu- 
rations NM, FM and AFM are given in table 5. 

The difference in values between the |J| of KNi[Cr(CN)g] 
and KV[Cr(CN) 6 ](40.5 cm" 1 for 296.5) is that the factor 
4/3S in the relation of J given by the Kahn's model J = 2k + 
4/3S [32 is equal to for the KNi[Cr(CN) 6 ] where it is neg- 
ative in the second one. k is the bielectronic exchange in- 
tegral (positive) between the two non orthogonalizcd mag- 
netic orbitals a and b, f3 is the corresponding monoelec- 
tronic resonance or transfer integral (negative) and S the 
monoelectronic overlap integral (positive) between a and 
b. It should be noted that 4/3S is more larger than 2k which 
gave a biggest value of |J| for the KV[Cr(CN) 6 ] complex. 

5. Conclusion 

The goal of this work was to carry out a theoretical 
study of the exchange coupling constants in two Prussian 
blues analogues as well as the magnetic order and chemical 
bonding. We chose a compound with a strong antiferro- 
magnetic coupling between the two TM, KV[Cr(CN)g] and 
with the same local spin, the other KNi[Cr(CN)g] with a 
ferromagnetic coupling. The sign of J, the values and the 
atomic spin densities obtained with the molecular calcula- 
tions are in good agreement with those obtained by other 
authors. The complex with t^-t^ pairs give a strong an- 
tifcrromagnctic interaction (case of Cr-V) while t^-t^c^ 
pairs yielding ferromagnetic coupling (case of Cr-Ni). This 
is an agreement with the conclusions given by Ruiz |13| . We 
note that the Tc values obtained from coupling constants 
in these two dinuclear models are significantly higher than 
the experimental ones. 

In the extended solid state approach, we have applied the 
all electrons ASW method to describe the electronic and 
magnetic properties as well as the chemical bonding of the 
two compounds. Our calculations yield stable and insu- 
lating magnetic ground states for both complexes. From 
the ECOV allowing for a qualitative chemical bonding 
analysis between atomic species, it has clearly been estab- 
lished that Cr must be bonded to C and that V to N. The 
Cr-C bond strength is found larger than the V-N one; a 
result confirmed by the vibrational spectra analysis. From 
the antiferromagnetic calculations, we proved that the 
KV[Cr(CN)@] complex is an insulating antiferromagnet 
while KNi[Cr(CN)g] is an insulating ferromagnet. 

The issue of the solid state calculations beside molecular 
ones is a major feature brought by this work by presenting 
the two approaches as complementary and bringing new re- 
sults such as casting a quantitative description for the en- 
ergy discrimination between different magnetic configura- 



tions and different chemical bonding possibilities through 
the isomerism. 
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Table 3 

Spin magnetic moments of each atom of the cell for the KV[Cr(CN)g] 
complex in the ferromagnetic (FM) and antiferromagnetic (AFM) 
configurations. ESI, ES2, ES3 are empty spheres entered within the 
ASA approximation (cf. text). 





FM 


AFM 


Atoms 


Magnetic moment value in 


Ms 


K 


0.005 


-0.001 


V 


2.409 


-2.273 


Cr 


2.681 


2.529 


C 


-0.006 


-0.076 


N 


0.008 


-0.045 


ESI 


0.004 


-0.041 


ES2 


0.028 


-0.047 


ES3 


0.042 


-0.066 


cell 


~ 6 


~ 0.0 



Table 1 

Antiferromagnetic and ferromagnetic exchange interactions between 
d orbitals of Cr and M ions. 



Cr CN 


M 


Type of interaction 


e g (empty) p CT 


e 9 (empty) 


AFM 


t2 g (half- filled) p„- 


t 2g (half-filled) 




e g (empty) p CT 


c g (partially filled) 


FM 


i,2g (half-filled) Pir 


t 2g (filled) 





Table 2 

Atomic spin densities on Cr, M for the two systems calculated with 
B3LYP functional and LanL2DZ basis set. 



System 


2S+1 Cr M 


Cr J "-V" 


7 3.10 2.57 


Cr J "-V" 


1 2.98-2.50 


Cr m -Nr" 


6 3.10 1.65 


Cr'"-Ni" 


2 3.08-1.65 



Table 4 

Spin magnetic moments of each atom of the cell or the KNi[Cr(CN)g] 
complex in the ferromagnetic (FM) and antiferromagnetic (AFM) 
configuartions. ESI, ES2, ES3 are empty spheres entered within the 
ASA approximation (cf. text). 





FM 


AFM 


Atoms 


Magnetic moment value 


in iiB 


K 


-0.000 


-0000 


Ni 


1.350 


-1.352 


Cr 


2.720 


2.715 


C 


-0.043 


-0.236 


N 


0.134 


-0.287 


ESI 


0.00 


-0.06 


ES2 


0.01 


-0.07 


ES3 


0.03 


-0.09 


cell 


~ 5 


~ 0.0 



Table 5 

Energy difference relative to the nonmagnetic state of the ferromag- 
netic and antiferromagnetic configurations for the two complexes, 
E NM = -85254.5512 eV for KV[Cr(CN) 6 ] and E NM = -100783.3848 
eV for KNi[Cr(CN) 6 ] 



A B (eV) 


V[Cr(CN) 6 ] 


Ni[Cr(CN) 6 ] 


NM 








FM 


-1.16 


-1.59 


AFM 


-1.31 


-1.51 




Figure 1. Sketch of the KM[Cr(CN) 6 ] structure. M stands for the 
other transition metal ion, V or Ni. 
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Infrared Spectrum 




SCO 1 COD 15D0 2OC0 250D 

Frequency (cm ) 

(b) 

Figure 2. KV[Cr(CN)e]: a) IR and b) Raman Spectrum calculated 
with B3LYP hybrid functional and LanL2DZ basis set. 
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Figure 4. KV[Cr(CN)6] in SP configurations: a) FM site projected 
density of states, b) AFM total DOS. 



Figure 3. KV[Cr(CN)6] in NM configuration: a) Site projected den- 
sity of states; b) Chemical bonding from ECOV criterion for two 
atom interactions: V-C, V-N, Cr-C and Cr-N. 
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